HEAT CONTENT ASYMPTOTICS WITH 
TRANSMITTAL AND TRANSMISSION BOUNDARY 

CONDITIONS 



^ ! P. GILKEY AND K. KIRSTEN 

o 

I Abstract. We study the heat content asymptotics on a Riemann- 

£^ ' ian manifold with smoooth boundary defined by Dirichlet, Neu- 

^ | mann, transmittal and transmission boundary conditions. 

Subject Classification: 58J50 



1. Introduction 



> 

(N . 

q \ Let M be a compact m dimensional Riemannian manifold with 

smooth boundary dM. Let D be an operator of Laplace type on a 
vector bundle V over M. Let B be a suitable local boundary condition 
and let D$ be the associated realization. Let <fi e C°°(V) describe the 
"*r^ \ initial temperature distribution. The subsequent temperature distri- 

Q-c bution u := e~ tDts (f) for t > is described by the equations: 



(1.1) {d t + D)u = 0, u(x; 0) = 0, and £w = 0. 
The specific heat p is a section to the dual bundle V*. Let 

f3(^p,D,B)(t) :=J M up 

j_i 1 be the total heat energy content. As t j 0, there is a complete asymp- 

totic expansion of the form 

(3(<P, Pl D, B){t) ~ E„> P> A B T'^ 

the /ieot content coefficients fl n ((f>, p, D, B) are locally computable. 

If Db is self-adjoint, then let {(pi, Aj} be a discrete spectral resolution. 
Let 7i(0) := J M 00« be the associated Fourier coefficients. Then: 

(1.2) P(<t>,p,D,B)(t) = Y,ie- tXi li{<t>)li{p)- 

It is convenient to introduce a formalism to consider both Dirichlet 
and Robin boundary conditions at the same time. Suppose given a 
decomposition dM = Cd U Cr of the boundary as the disjoint union of 
two closed (possibly empty) sets. Let S be an auxiliary endomorphism 
of V\c R and let ;m be the covariant derivative of </> with respect to 

l 
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the inward unit normal, where we use the natural connection which is 
induced on V by D - see Section for details. We define 

B D r = B d ®B r where B D (j) := (j)\c D and B R (f) := (<p. m + S(j))\c R 

are the pure Dirichlet and Robin operators respectively. In Section |2|, 
we review previous results for the boundary conditions Bdr- 

Transmittal and transfer boundary conditions will form the primary 
focus of this paper. Let (M±,g±) be smooth compact m dimensional 
Riemannian manifolds. We assume that E = dM + = dM_ is a smooth 
m — 1 dimensional manifold and that the induced metrics agree, i.e. 
<7+|s = (7— |s- Let D± be operators of Laplace type on vector bundles 
V± over M± . Let v± be the inward unit normals of E C M± ; note that 
u + = —V-. Let (ft := (0 + , 0_) and p := (p + , p_). 

Suppose that V+|s = F-ls and that there is given an auxiliary en- 
domorphism U of Vs := V±|s serving as an impedance matching term. 
Let V be the natural connections defined by the operators D±. Let 

Bi4>: = {0+|s-0-|e} 

(1.3) © {(v+ + + )| s + (v;_0_)| s -f/0 + | s }. 



Equivalently, satisfies the boundary conditions given in display ( |1.3|) 
if and only if extends continuously across the interface E and if the 
normal derivatives match, modulo the impedance matching term U. In 
Section || we determine the invariants f3 n for n < 3 for these bound- 
ary conditions, see Theorems [3.i| and |3.2| for details. The transmittal 
boundary operator B\ = B\(U) is of relevance in the presence of distri- 
butional sources [K| |13| as they have been considered, e.g., in the 
brane world scenario. 

We shall also be studying boundary conditions which are defined by 
the boundary operator i3 2 = ^(5): 



(1.4) B 2 
where 



v+ + S ++ \(<j )+ 
S-+ V- J [ 0_ 



If 5 + _ = S'_ + = 0, then equation ( |1.4| ) decouples to define Robin 
boundary conditions. Note that we do not assume given an identifi- 
cation of V+|s with V-.\y,; in particular, we can consider the situation 
when dim V + ^ dim V- . In Section f| we determine the heat content 
invariants j3 n for n < 3 for the heat transfer boundary conditions B 2 , 
see Theorem 14.31. 
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The boundary conditions defined by equations ( |1.3| ) and ( |1.4j ) can be 
thought of as living on the singular manifold M := M + Ue M_. Both 
boundary conditions are relevant to heat transfer problems between 
two media of different conductivities. Which boundary condition is 
to be applied depends on the details of the surface of separation £ 
between M + and M_. Let K + and K_ be the thermal conductivities 
of M + and M_. The flux of heat is continuous over the interface S, 

(i.5) (ir+v+^ + ^-V- =0. 

If the contact between the two media M + and M_ is very intimate, in 
addition one assumes 

(1-6) 0+| E = 0_| s , 

and boundary conditions of the type (|1.3|) are found. Otherwise, e.g., 
for surfaces pressed lightly together, in a linear approximation the flux 
of heat between M + and M_ is proportional to their temperature dif- 
ference. In this case, equation (|1.5| ) has to be augmented by 

(1-7) (K + V+ + 0+)| S =i?(0 + -0-)| s 

where H is referred to as the surface conductivity. The boundary 
conditions ( |1 . 5[ ) and fll.Tj ) can be combined into the form of equation 
(|1.4|); see, for example, the discussion in ||. 

2. DlRICHLET AND ROBIN BOUNDARY CONDITIONS 

We begin by reviewing some of the basic invariance theory of opera- 
tors of Laplace type. Let (M, g) be a compact Riemannian manifold of 
dimension m. We suppose the boundary £ of M is smooth. We adopt 
the Einstein convention and sum over repeated indices. Let 

D = -(g^d fl d u + A< i d tl + B) 

be an operator of Laplace type on C°°(V). The operator D determines 
a natural connection V and a natural endomorphism E such that we 
may express D invariantly in the form: 

D = -{Tr(V 2 ) + E}; 

see || for details. Let T be the Christoffel symbols of the metric. We 
may express the connection 1 form u of V and the endomorphism E: 

u 5 = \gus{A v + g^T,^) and 
(2.1) E = B — g^id^ + uy^ - u a V u ^). 

Note that the connection defined by the dual operator D on the dual 
bundle V* is the associated dual connection with connection 1 form 
given by — uj*; furthermore the associated endomorphism is E*. 
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We shall let Roman indices a, b, etc. range from 1 to m — 1 and 
index a local coordinate frame for the tangent bundle of the boundary. 
Let e m be the inward unit normal and let indices i, j, etc. range 
from 1 to m and index this augmented frame for TM. Let L ab be the 
second fundamental form, let Rijki be the Riemann curvature tensor 
with the sign convention that R1221 = +1 for the unit sphere in R 3 . 
Let Q be the curvature of the induced connection on V. Let ':' and ';' 
denote multiple covariant differentiation with respect to the Levi-Civita 
connection of the boundary and of the interior, respectively; these two 
connections differ by the second fundamental form. 

The invariants j3 n may be decomposed as sums j3 n = + (3^ of 
locally computable invariants given by integrals over the interior and 
over the boundary. Let D and B be the dual operators on C°°(V*). 
The interior invariants are independent of the boundary condition and 
vanish if n is odd. For n < 3, we have: 

$*{p, <f), D, B) = j M 4> ■ p, <P, D, B) = 0, 

$*{p, 0, D,B) = - J M Dc/> ■ p, /3f 0, D, B) = 0. 

The heat content asymptotics (3 n defined by the Dirichlet and Robin 
boundary operator Bdr have been studied previously [J], |2], || J|, II], [12 



[L4], [L^, [HJ. There are also results available in the singular setting, see 
for example || We summarize the results for /3 , /3i, /3 2 , and /3 3 : 

Theorem 2.1. 

1. f3 ((f),p,D,B DR ) = f M <f>p. 

2. fr(<l>,p,D i B DR ) = -%f CD <l>p. 

3. (3 2 ((f), p, D, B DR ) = - J M D(p- p + j CD {\L aa cj)p - 0p ;m } 

+ J Cr Br<P ■ p. 

4. p, D, B DR ) = -^j Co {-lD4> ■ p - l<j>Dp + l^ap-.a 

{—\E + j^L aa L bb — ^L ab L ab + ^R a mam)4>P} +3^ / C R ^ R( ^ ' BrP- 

3. The boundary operator B\ 

We postpone for the moment the discussion of f3 3 . Using the chiral 
symmetry and the homogeneity of the invariants, we see: 

Theorem 3.1. There exist universal constants so 

1. p {(f), p, D, B{) = j M+ <p + p + + j M _ . 

2. /3i(c/>, p, D, B x ) = / E {ai(0+p+ + (j>-p-) + a 2 ((f> + p- + <f>-P+)}- 

3. (3 2 {(p, P, D, B{) = - j M+ D + <p + ■ p + - f M _ £L0_ ■ p_ 

+ k { a z(<t>+P+ L ia + <P-P- L bb) + a i(.<P+P+ L aa + <P-P- L ta) 
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+a 5 (0+p_L+ + <j>-p+L- a ) + a 6 (0+p_L- + (p-p+L+ a ) 
+a 7 ((p +;u+ p + + 4>-- u _p-) + a 8 ((/> +; „ + p_ + (f)- tU _p + ) 

+a 9 (4> + p + - y+ + 4>-P-;vS) + aio(0+P-;i._ + 

+an(0 + p+ + (p-P-)U + ai 2 (0+p- + 0_p + ){7}. 
4. VFe /iawe: 

a i = — a 2 = o 3 = |, a 4 = |, a 5 = — ~, a 6 = — ~, 

07 = |, a§ = |, ag = — |, aio = |, an = — j, ^12 — — j- 

There are a number of functorial properties that these invariants 
satisfy. Suppose that the bundles V± are equipped with Hermitian 
inner products, that the operators D± are formally self-adjoint, and 
that U is self-adjoint. We then have that D is self-adjoint, see 
(equation (17)) for details. We may therefore apply the relations of 
display ( |1.2| ) to see: 

(3 n ( ( p,p,D,B 1 )=(3 n (p, 0,A#i)- 

More generally, if D is the formal adjoint of D on C°°(V*) and if B\ 
are the dual boundary conditions, then we have 

(3.1) (3 n ((P,p,D,B 1 )=/3 n (p,(P,D,B 1 ). 

The expression — f M Dtp ■ p is not symmetric in tp and p. We use 
equation (|3.1|) and integrate by parts to see: 

(3.2) a 5 = a 6 , a 7 - a 9 = 1, a 8 = a w . 

Doubling the manifold yields additional information. Let M be a 
smooth Riemannian manifold with smooth boundary £ and let D be 
a self-adjoint operator of Laplace type over M . Let {<pD,i,^D,i} and 
A^j} be the discrete spectral resolutions for D with Dirichlet 
(D) and Robin (R) boundary conditions over M Q . Let M ± := M 
define the double. Extend the <pD,i to be odd and the 4>r,% to be even: 

<pD,i(x±) = ±^tp D ,i( x ) and <t>R,it X ±) = ^4>R,ii X )- 

Set U = —2S. It was shown in [|K| that B\<pD,i = and B\4>R,i = 0. 
Furthermore, {(pD,i, <pR,i} is a complete orthonormal basis for L 2 (V) 
which defines the spectral resolution of D := (Dq , Dq). Decompose 
4> = <Po + <Pe and p — p + p e as the sum of even and odd functions and 
let tp , (p e , p , and p e be the restrictions to M = M + . We then have 

4>o = Y,ilD,i(4>o)<pD,i, <Po = V2Y,ilD,i{<t ) o)<pD,u 

Pp = E 4 lD,i(po)(pD,i, Po = V^Ei lD,i{Po)<t>D,h 

= Ei7R,i(0e)^ii,i, 0e = Ej 7fl,i (0e) 0H,i , 

Pe = Ei lR,i(P~e)4>R,i, Pe = V 7 ^ Ei lR,i(pe)<pR,i- 
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Consequently by equation ( |1.2| ), 

(3(<t>, p, D, BJit) = 2(3(4> , p , D , Bo) (t) + 2/3(} e , p e , D , Br) (t) 

(3.3) /3 n (0, p, D, Bt) = 2f3 n (4> , p , D , B D ) + 2/3 n (0 e , p e , D , B R ). 

This relation continues to hold even if Dq is not self-adjoint. Thus 

2ai + 2a 2 = 0, 2d\ — 2a 2 = — 

2a 3 + 2a 4 + 2a 5 + 2a 6 = 0, 2a3 + 2a 4 — 2as — 2a6 = 1, 

(3.4) 2a 7 + 2a 8 = 2, 2a 7 - 2a 8 = 0, 
2a 9 + 2a 10 = 0, 2ag — 2aio = — 2, 
— 4an — 4ai2 = 2, 2an — 2a\ 2 = 0. 

Take arbitrary metrics on M± and let D± be the scalar Laplacian. 
Take <p = 1 and U = 0. Then £>0 = and #i0 = so e- tD ^(j) = <p. 
Thus /3 n (l, p, D, B\) = for n > 1. Take p_ = 0. The terms p + , p+L+ a , 
p+L~ a , p+- u+ can then be specified arbitrarily. This yields: 

(3.5) a\ + a,2 = 0, CL3 + as = 0, (Z4 + 05 = 0, ag + aio = 0. 

This allows for the determination of the multipliers a±, 0,12. How- 
ever, in order to provide further checks and because it will be useful 
later, we give one final property. Let N± := [0, 1] be the interval. Let 
M± := [0, 1] x S* 1 be the cylinder with the metrics 

ds 2 = dr 2 + e 2f ± (r) dd 2 

where the real functions f± vanish on <9{[0, 1]}. Let f± yT := d r f±. Let 

D ±)N := —(d* + f±, r d r ) on N± and 

D ± ,m = -(d 2 r + f±,rd r + e~ 2f± d 2 e ) on Mfc. 

Then D±^m is the scalar Laplacian on M±. The second fundamental 
form vanishes on N± while = —f± r is the second fundamental form 
on M±. The connection forms defined by these two operators differ: 

= \f r on V, = -\f r on V*, 
= on V, = on V*. 

To compensate for this difference, we let 

U N = i(/ +ir + /_ r ) on ON and U M = on <9M. 
The volume forms also differ: 

dvol N = dr on N± and dvol M = e f± drd9 on M±. 
We let 4>± and p± be constants. We then have: 

e * Db i cj) = e * Db i so 

(3.6) f3 n ((f>,p,D M ,BM) = 2Tr{3 n ((t),efp,D N ,B?). 
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We take f + = f and /_ = 0. On the cylinder, the only invariant 
that plays a role in the computation of /?2(0, P, D, B\) is L+ a = —f r . 
On the interval, the only invariants that play a role are U = \ f ri the 
connection 1 form u) r = |/ r on V, the dual connection one form —\f r 
on V*, and the endomorphism E = —\fr~\frr\ the interior invariants 
vanish as 

D M ( ( j ) ) = 0, D N (cf>) = 0, D M (p) = 0, D N (efp) = 0. 
Note that on dN, (e? p + ) ]U+ = ^f r p+- Thus equation ( |3.6|) implies: 
f r (-a 3 (p + p + - a^-P- - a 5 <p + p^ - a 6 0-P+) 
= \fr{ai4>+p + + a 8 + p_ + a 9 cf) + p + + a lo 0_p + 
+an(0 + p + + 0_p_) + ai 2 (0+p- + 4>-P+)} 

and consequently 

/o y\ — 2a 3 = 07 + 09 + an, — 2a 4 = an, 

' — 2a 5 = a 8 + ai2, — 2a 6 = a 10 + ai 2 - 

We solve the relations of displays fl3.2|) , (|3.4j) , fl3.5|) , and (|3.7|) to com- 
plete the determination of /3q, and /?2 in this setting by determining 
the unknown coefficients to complete the proof of Theorem |3.1| . □ 

Let u a : = V+ — on V and £o a = — u* a on V*; this is a chiral tensor 
that changes sign if we interchange the roles of ± or of V and V*. We 
determine /3 3 in this setting: 

Theorem 3.2. 

1. There exist universal constants so /3 3 (0, p, D, B%) 

= 67^ J s { a 2o{D+<f>+ ■ P+ + 0+ ■ D + p + + LL0_ ■ p_ + 0_ ■ D-p-) 

+a 2 i (£>+</>+ • p- + 0+ ■ D-p- + LL0_ ■ p + + 0_ • 

+a 22 (^ a V+0+ ■ p+ - u; a V^0_ ■ p- - u a (p + ■ Y+p+ + W a 0_ ■ V« P-) 

+a 23 (^ a V+0+ -J- - LO a V~(j)- ■ p+ + U a (j) + ■ V~p- - W a 0_ ■ V+p+) 

+a 24 (V+ 0+ ■ V+p+ + V;_0_ • V-_ P -) 
+a 25 (V+ 0+ • V;_p_ + V;_0_ ■ V+ p+) 
+a 26 (V+0+ ■ V+p+ + V"0_ • V"p_) 
+a 27 (V+0+ • V-p_ + V~0_ • V+p+) 
+a 28 [/(<% + (0+p + ) + <%_(0_p_)) 

+a 29 [/(V;_0_ ■ P+ + 0- ■ V+p+ + V+0+ ■ p_ + 0+ • V- p_) 

+«3o(^aa^+(0+P+) + L aa d v-{<l>-P-)) 

+azi{L- a d u+ ((j) + p + ) + L+ d v _{(j)-p-)) 

+a 32 (L+ (V+ 0+ • p- + 0- V+ p+) + L- a (V;_0_ ■ p+ + 0+V;_p_)) 
+a 33 (L- a (V+ 0+ • p_ + 0-V+ p + ) + L+ (V- 0_ ■ p+ + 0+V;_p_)) 
+a 34 tj a ^ a (0 + p + + 0-P-) + a 35 w a a; a (0 + p_ + 0_p+) 
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+«38 (£aaW+P+ + LtaLtb^-P-) + a 39 L ia L bb( < l ) + P- + 

+a 40 (L+ + L~ a L bb )((j)+p-. + (f)-p+) 

+a 41 (L+ b L+ b (f) + p + + L ab L ab 4>_p„) + a 42 L+ b L ah (<j)+p+ + <j>-p-) 

+a A3 (L ab L ab cf> + p+ + L+L+ b <f)„p_) + a 44 L+ b L ab ((f) + p_ + 

+a 45 (L+L+ + L ab L ab )(c/) + p_ + <p-p + ) _ 

+a m U{L+ a (j) + p + + L aa 0_p_) + a 47 [/(L aa 0+p + + ^aa^-P-) 

+a 48 C/(L+ + L~ )(0+p- + <p_p + ) 

+a A9 U 2 ((f) + p + + 4>-P-) + a 5Q U 2 ((f> + p- + 4>-P+) 

+a 51 (E + cf) + p + + £_(/>_p_) + a 52 (E-<f>+P+ + 

+a 53 (E + + Ej){<j> +P _ + 0_p + ) + a 5A (R± ji (j)+p + + R'^-pJ) 

+att(.Rjjji<P+P+ + Rtjji<f>-P~) + a w{ R tjji + R ijji)(<P+P- + <P-P+) 

+ a 57 (Rtmma ( fi+P+ + ^ amma 

+^58 (Ramma C f ) +P+ + Rtmma^-P-) 

+Q.59 (R amma + R a mma ){4>+p- + 4>-p+)- 



2. a 20 = 


4, 


«21 = 


-4, 


«22 — 


-1, 


«23 = 


-1, 


&24 = 


4, 


^25 = 


4. 


026 = 


-2, 


a 27 = 


2. 


028 = 


-2, 


«29 = 


-2, 


«30 = 


-1, 


^31 = 




«32 = 


1. 


«33 = 


-1, 


«34 = 


1, 


«35 = 


0. 


036 = 


0. 


«37 = 


1 

2' 


a 3 8 = 


0. 


a 39 = 


1 

2 • 


a 4 o = 


0. 


a 4 i = 


1 

2 1 


a 4 2 = 


0. 


&43 — 


1 

2 1 


CZ 44 = 


0. 


a 4 5 = 


1 

2' 


a 46 = 


1, 


a 4 7 = 


-1, 


a 48 = 


0, 


a 49 = 


1. 


&50 = 


1. 


051 = 


1. 


^52 = 


1. 


053 = 


-1, 


^54 = 


0, 


&55 = 


0, 


056 = 


0. 


057 = 


1 

2 - 


«58 = 


1 

2 1 


a 59 — 


1 

2' 



We use the relations of equation ( |3.1| ) to simplify the format at the 
outset and derive (1). We shall use the functorial properties involved 
to determine the unknown coefficients and prove (2). 

We apply Theorem |2J] and equation (|3.3|) to see: 



2a 2 o + 2a 2i = 0, 
2a 24 + 2a 25 = 16, 
2a 26 + 2a 27 = 0, 
-4a 28 - 4a 29 = 16, 
a 3 o + a 3 i + a 32 + a 33 = 0, 
a 3 6 + a 3 7 + a 3 s + a 3 9 + 2a 40 
a 4 i + a 42 + a 4 3 + a 44 + 2a 4 5 
2a 46 + 2a 47 + 4a 48 = 0, 
8a 49 + 8a 50 = 16, 
2a 5 i + 2a 52 + 4a 53 = 0, 
2a 54 + 2a 55 + 4a 56 = 0, 
2a 57 + 2a 58 + 4a 59 = 0, 



0, 
0, 



2a 20 — 2a 2 i = 16, 

2a 24 - 2a 25 = 0, 

2a 2 6 — 2a 27 = —8, 

2a 28 - 2a 29 = 0, 

2a 30 + 2a 3 i - 2a 32 - 2a 33 = 0, 

a 36 + &37 + a 38 — a 39 — 2a 4 o = — 1, 

a 4 i + a 42 + a 43 — a 44 — 2a 4 5 = 2, 



2a 4 6 


+ 2a 47 


- 4a 48 


= 


2a 49 


— 2a5o 


= 0, 




2a 5 i 


+ 2a 52 


- 4a 53 


= 8 


2^54 


+ 2a 55 


- 4a 56 


= 


2a 57 


+ 2a 58 


- 4a 59 


= 4 
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Take arbitrary metrics on M± and let D± be the scalar Laplacian. 
Take = 1 and U = 0. Then (p satisfies transmittal boundary condi- 
tions. Thus /3 n (l, p, D, Bi) — for n > 1. Take p_ = 0. This yields: 

^30 + a 3 2 = 0, a 3i + a 33 = 0, 

a 36 + a 4o — 0, a 3 g + = 0, a 37 + a 3 g = 0, 

a4i + 045 = 0, a4 3 + 045 = 0, 042 + 044 = 0, 

a 5 i + «53 = 0, a 52 + a 53 = 0, a 54 + a 56 = 0, 

«55 + ^56 = 0, a 57 + a 59 = 0, a 58 + a 59 = 0. 

Let D±{e) = D± - e. Then D±(e) = D± - e and E±(e) = E± + e. 
As e ~ tDB ^ = e t£ e- tD ^, P(<f>,p,D(e),Bi)(t) = e te (3((j), p, D, B\){t) and 
hence d £ j3 n \ e= o = (3 n - 2 . Thus studying the coefficients of the terms 
{4> + p + , 4>+P-} leads to the relations: 

g^{— 2a 2 o + a 5 i + a 52 } = a x = 

6^{-2a 21 + 2a 53 } = a 2 = 

We use separation of variables to generate additional relationships 
among the coefficients. First, we study flat metrics. Let (r, 9) be the 
usual parameters on M± := [0, 1] x S 1 . Let 

D± := -(d 2 + d 2 + 2e ± d e ) 

where e± = e±{6). Let N± := [0,1] and let D± := —d 2 . Let <p± and 
p± be constant. Let Uq be constant. Separation of variables and an 
application of equation ( |3.3| ) and Theorem [^T] then yields: 

e - tDB i(J) = e~ tD± (j) so 

M, p, D, B l ) M = 27r/5 3 (0, p, D, B X ) N 

= 47r/? 3 (0e, Pe, 5, B R ) N+ + 47T/3 3 (0 O , p , D, B d )n+ = 0. 

As Ve0± = £±<j>±, Vgp± = — £±p±, and i?± = — e± — dgE±, we have: 

= a 22 {s + - e^){2s + (j) + p + - 2e„0_p_) 
+a 23 (e+ - £_) 2 (0+p_ + 

+a 26 (-£+0+P+ - £ 2 + a 27 £+£„(-0+p_ - 0-p+) 

+a 34 (e + - £_) 2 + p + + + a 35 (e + - £_) 2 (0 + p_ + </>_p+) 

+a 51 (-e 2 + (f) + p + - £ 2 4>~p-) + a 52 (-e 2 _(j) + p + - e\4>-p-) 
+a 53 (-4 - e 2 _ )(0+p_ + 

The terms {e 2 + p + , £ 2 + p + , e + e_0 + p + , e 2 + p_, £ + e_^ + p_} are then 
studied to conclude that: 

= 2a 22 — a 2 Q + a 3 4 — 051, = a 3 4 — a^ 2 , 

= -2a 22 - 2a 34 , = a 23 + a 35 - a 53 , 

= -2a 23 - a 27 - 2a 35 . 
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We can also get information from the divergence terms. We now let 
cf) = (1, 1) and p = (p + (9), 0). We work modulo 0(e 2 ) and use the fact 
that a 2 o + 0,21 = to show: 

= is {{dep+){-a>22(e+ ~ £-) - a 2 3(s+ - £-) + a 26 £+ + ^27^-} 

+p + {-a 51 d e e + - a b2 d e e- - a 53 d e (e + + £-)}} 
-a 22 - a 23 + a 26 + a 5 i + a 53 = 0, 
a 22 + a 2 3 + a.27 + a 52 + «53 = 0. 
Next, let N± := [0, 1] and M± := [0,1] x S 1 x S 1 have the metrics: 



ds 2 = dr 2 + e 2h,±(r) d0 2 + e 2f 2 Mr) M 



where f^± vanishes on the boundary. Let fi,±, r '■= d r fa,± an d let 

d£ :=-{fl? + EiA±A}. 

£>£ := -{d 2 r + Y,iUi,±A + e-VwflJ)}. 

We set t/jvf = i/o constant. As the connection forms defined by D N and 
D M are different, we set Un = Uo + (/*,+,»• + fi,-,r)- Let 0± and 
p± be constant. The argument used to establish equation (|3.6| ) then 
generalizes immediately to yield: 

(3.8) M,P,D M ,B 1M ) = (27i) 2 p 3 (<p,e^ f 'p,D N ,B 1 , N ). 



We compute 


on <9M: 






r ± 

mab 


— 2 c 'mi/a6 ~~ 






r ± 

abm 


- -1/3 o ± 


^ "'^ fa,±,r$abi 




r ± 

amb 


- -T^ - 

abm 


& ^ ' f a,±,r&abi 




p± b 

am 


— y 1 amb ~ 


fa,±,r"abi 




T± 
^ab 


= (V a <9 b , <9 m 


\ —V ± — - P 2 f*,± f 

) ~ 1 abm ~ c V a 




D± 

-^amma 


= f(V ± V ± ■ 


-V±V±),d m; d a ) = 


_ (V m r am db,d m ) 



— (V m / a ,±,r<9 a , l9 a ) — fa,±,rr fa,±,V 

(We do not compute iE-^ as we have shown a 54 = a 55 = a 56 = so 



these terms do not appear). Let f± = ^ /^-t. We compute on ON: 

Un = Uq + 2(/+,r + /-,r), -E'at = — ■ 2 ~f±,rr ~ lf±,r i 

Vtp± = \f±4>±, V±(e'±p) = ±/±p±, 

D£<f> ± = 0, DUe f± P±} = 0. 



We use equation ( |3.8j ) to derive the following equations from the coef- 
ficients of the indicated monomials: 
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We now let p = p(r) . Equation ( |3.8| ) continues to hold and yields: 



We combine the relations given above to complete the determination 
of p 3 by determining the constants a 20 — a 59 and complete the proof of 
Theorem 13.21. □ 



Again, we begin our discussion by studying p , Pi, and Note 
that terms such as cannot occur since we do not assume an 

identification of V+ |s with V- |s- 

Lemma 4.1. There exist universal constants so 



1. p ((f), p, D, Bi) = j M+ <f) + p + + J M _ (f)-p- . 

2. ft(0,p,£,B 2 ) = / E & 1 (0+P+ + 0-P-)- 

3. /3 2 (0, p, D, B 2 ) = - J M+ D + <p ■ p + - f M _ ■ p_ 

+ fxih^+p+L+a + (j>-P-L bb ) + b 3 ((f) + p + L- a + (f>-P-L+ a ) 

+h(<p +]U+ p+ + (j)-,v_p-) + h((p + p + . u+ + 4>-p-. u _) 

+b 6 (S ++ (f> + ■ p + + • p_) + 6 7 (5 + _0_ • p + + • p_)}. 



~ a 31 — |&28 

— ^33 — 2 fl 25 + 2 fl 29 



(fl,-r(i>+d„ + p + ) 
Ul,+,r<P+d u _p-) 



4. The boundary condition B 2 
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Taking S' + _ = and S-+ = forces the boundary conditions given 
in equation (|1.4j) to decouple and defines Robin boundary conditions 
on M + and on M_ separately. We use Theorem [O] to see: 

(4.1) b 1 = 0, b 2 = 0, b 3 = 0, 6 4 = 1, h = 0, 6e = l. 

We use an argument similar to that used to establish display ( |3.5|) to 
determine b-j. Let D± be the scalar Laplacians on manifolds M±. Let 

4> + = 4>- = 1, let S ++ = S = 1, and let S+- = S_ + = —1. Then 

B 2 <p = and D<p = so (3 n = for n > 0. Thus: 

(4.2) b 6 -b 7 = 0. 

In view of the remarks noted above, we see that: 
Lemma 4.2. There exist universal constants so 
P 3 ^ 1 p,D,B 2 ) = ^J^B 2( f ) -B 2 p 
+& lo (5+_5_ + + • p+ + 0_ • p_) 

+b n (S—S-+<f>+ ■ p- + S ++ S + _</>_ • p+) 
+b 12 (S_ + S ++ (f) + ■ p_ + S , + - 5 '— 0- • p + ) 
+b 13 (S- + (f} + - u+ ■ p- + S+-<t>-. v _ ■ p+) 
+b 14 (S^ + (f) + ■ p„. v _ + S+-(f>- ■ p+-u+) 

+ bl5( L ta S - + ( l ) + ■ P- + L aa S +~ ( t ) ^ ' P+) 

+6 16 (X- ■ p_ + L+ S+_0_ • p+). 

Let D-t be the scalar Laplacians on manifolds M±. Let + = 1, 

4>- = 1, let S' ++ = a, S^- = —a, 5 = 6, and >S_ + = —6. Then 

satisfies transmittal boundary conditions and is harmonic so f3 n = 
for n > 0. Consequently taking p_ = yields the equations: 

(4.3) 610 - 612 = 0, b n = 614 = 615 = 616 = 0. 

We work with m — 1 and D± = —d%. Suppose that <fi± = a±x + b± is 
such that 4> satisfies B 2 <p = 0, i.e. 

(4.4) ea + + S ++ b + + S+-b_ = and ea_ + + SL+&+ = 

where e(0) = +1 and e(l) = —1. We choose 5* and 6* arbitrarily and 
use equation ( |4.4|) to determine a*. Let p_ = 0. Since /3 n = for n > 0, 

b w S + -S-+b + + b 12 S + -S — + bi 3 £S + -a_ = 0, 
(610 - b 1 3)S+S- + b + + (612 - f>i 3 )5' + _S'__b_ = so 

(4.5) 610 = 613, and b 12 = b 13 . 
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We double the manifold to complete our determination. Suppose 
given an operator Do of Laplace type on a manifold Mo with bound- 
ary S. Let an initial condition <f>Q be given and let uq solve equation 
fll.lj) with the boundary operator Br. Let M± := M and D± := D . 
Then u± := u and (j>± := 4> solves equation ( |1.1| ) with the bound- 
ary operator B 2 defined by equation ( |1.4| ) with S ++ = S = and 

S + _ = <5_ + = S. Thus 

(3 n ((j)o,p + + p^,D ,B R ) = p n (^,P,D,B 2 ). 

We may now conclude b 10 = b 13 = 0. This proves 

Theorem 4.3. 

1. o ((t>,P,D,B 2 ) = f M (t)p. 

2. p 1 (<t>,p,D,B 2 ) = 0. 

3. /3 2 (0, p, A ft) = -J M D(j>-p + fx B 2 <f> ■ p. 

4. (3 3 (ct>,p,D,B 2 ) = 1 ^f^B 2( l>-B 2 p. 
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